Negative conductances of Josephson junctions: Voltage fluctuations and energetics 
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We study a resistively and capacitively shunted Josephson junction, which is driven by a com- 
bination of time-periodic and constant currents. Our investigations concern three main problems: 
(A) The voltage fluctuations across the junction; (B) The quality of transport expressed in terms of 
the Peclet number; (C) The efficiency of energy transduction from external currents. These issues 
are discussed in different parameter regimes that lead to: (i) absolute negative conductance; (ii) 
negative differential conductance, and (iii) normal, Ohmic-like conductance. Conditions for optimal 
operation of the system are studied. 
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I. INTRODUCTION 

Transport processes in periodic systems play an im- 
portant role in a great variety of everyday life phenom- 
ena. Two prominent examples are the electric transport 
in metals providing a prerequisite of modern civilization, 
and the movement of so-called molecular motors (like ki- 
nesin and dynein) along microtubules in biological cells 
which are of crucial relevance for the functioning of any 
higher living organism. Josephson junctions belong to 
the same class of systems being characterized by a spa- 
tially periodic structure. In the limiting case of small tun- 
nel contacts the mathematical description of a Josephson 
junction is identical to that of a Brownian particle mov- 
ing in a periodic potential. Such models have also fre- 
quently been employed under nonequilibrium conditions 
to describe Brownian ratchets and molecular motors, see 
Refs [1, 2] and references therein. Of particular impor- 
tance for technological applications are 'rocked' thermal 
Brownian motors operating either in overdamped or un- 
derdamped regimes, e.g. see in Ref. Q. The majority of 
papers on transport in periodic systems are focused on 
the asymptotic long time behavior of averaged quantities 
such as the mean velocity of a molecular motor, or the 
mean voltage drop in a Josephson contact [3] . The main 
emphasis of these works lies on formulating and exploring 
conditions that are necessary for the generation and con- 
trol of transport, its direction, and magnitude as well as 
its dependence on system parameters like temperature 
and external load. Apart from these well investigated 
questions other important features concerning the qual- 
ity of transport though have remained unanswered to a 
large extent. The key to these problems lies in the inves- 
tigation of the fluctuations about the average asymptotic 
behavior Q. 
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In the present paper we continue our previous studies 
on anomalous electric transport in driven, resistively and 
capacitively shunted Josephson junction devices 
These investigations were focused on the current-voltage 
characteristics, in particular on negative conductances. 
In contrast, in the present paper we investigate the fluc- 
tuations of voltage, the diffusion processes of the Cooper 
pair phase difference across a Josephson junction as well 
as the energetic performance of such a device. 

The paper is organized as follows. In the next section, 
we briefly describe the Stewart-McCumber model for the 
dynamics of the voltage across a junction. In Section 
3, we study voltage fluctuations, phase difference diffu- 
sion, and the efficiency of the device. Conclusions are 
contained in Section 4. 



II. MODEL OF RESISTIVELY AND 
CAPACITIVELY SHUNTED JOSEPHSON 
JUNCTION 

The Stewart-McCumber model describes the semi- 
classical regime of a small (but not ultra small) Josephson 
junction for which a spatial dependence of characteristics 
can be neglected. The model contains three additive cur- 
rent contributions: a Cooper pair tunnel current charac- 
terized by the critical current /q, a normal (Ohmic) cur- 
rent characterized by the normal state resistance R and 
a displacement current due to the capacitance C of the 
junction. Thermal fluctuations of the current are taken 
into account according to the fluctuation-dissipation the- 
orem and satisfy the Nyquist formula associated with the 
resistance R. The quasi-classical dynamics of the phase 
difference (p = (j){t) between the macroscopic wave func- 
tions of the Cooper pairs on both sides of the junction is 
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described by the following equation 0, [T^l , 

+Ia coa{nt + tpo) + , 



(1) 



where the dot denotes the differentiation with respect to 
time, Id and la are the amplitudes of the applied direct 
(dc) and alternating (ac) currents, respectively, O is the 
angular frequency and tpo defines the initial phase value 
of the ac-driving. Thermal equilibrium fluctuations are 
modeled by zero-mean Gaussian white noise ^(t) with the 
correlation function < ^(i)^(s) {2kBT/R) 5{t - s), 
where Ub is the Boltzmann constant and T is tempera- 
ture of the system. 

The limitations of the Stewart-McCumber model and 
its range of validity are discussed e.g. in Sec. 2.5 and 
2.6 of Ref. There are various other physical sys- 

tems that are described by Eq. ([1]). A typical example 
is a Brownian particle moving in the spatially periodic 
potential U{x) — U{x + L) = — cos{x) of period L — 2n, 
driven by a time-periodic force and a constant force Q. 
In this case, the variable corresponds to the spatial co- 
ordinate X of the Brownian particle and ac and dc play 
the role of periodic driving and a static tilt force, respec- 
tively, acting on the particle. Other specific systems are: 
a pendulum with an applied torque Q , rotating dipoles 
in external fields [ll|, Il2l . superionic conductors jl3l | and 
charge density waves |14| . 

It is convenient to transform Eq. ([T]) to a dimension- 
less form. We rescale the time t' = Upt, where ujp = 
{l/h)^/8EjEc is the Josephson plasma frequency ex- 
pressed by the Josephson coupling energy Ej — (ft/2e)/o 



and the charging 
takes the form [9 



energy Ec 
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V2C. Then Eq. (P) 



+ 7^ + sin((/)) =iQ + ii cos(riit' -I- ipo) 
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The dimensionless damping constant 7 = l/ujpRC is 
given by the ratio of two characteristic times: tq — l/ojp 
and the relaxation time Tr = RC . This damping con- 
stant 7 measures the strength of dissipation. The ac 
amplitude and angular frequency are ii = la/Io and 
VLi = U,tq — ^/ujp, respectively. The rescaled dc strength 
reads ia — Id/Io- The rescaled zero-mean Gaussian 
white noise T{t') possesses the auto-correlation func- 
tion {r{t')r(u)) — S{t' — u), and the noise intensity 
D — ksT/Ej is given as the ratio of two energies, the 
thermal energy and the Josephson coupling energy (cor- 
responding to the barrier height). 

Because Eq. ^ is equivalent to a set of three au- 
tonomous first order ordinary differential equations, the 
phase space of ^ is three-dimensional. For vanishing 
diffusion constant, D = Q, the system becomes deter- 
ministic. The resulting deterministic nonlinear dynam- 
ics [D = 0) exhibits a very rich behavior ranging from 
periodic to quasi-periodic and chaotic solutions in the 



asymptotic long time limit. Moreover, there are regions 
in parameter space where several attractors coexist. In 
the presence of small noise these attractors still domi- 
nate the dynamics in the sense that most of the time the 
trajectory stays close to one of these attractors. Only 
rarely, transitions between the attractors take place. So, 
the locally stable states of the noiseless dynamics become 
metastable states in the presence of weak noise. Apart 
from that, the presence of noise may also let the system 
come close to deterministic unstable orbits which it may 
follow for quite some time. 

Strictly speaking, the deterministic regime D = is 
only reached in the limit of zero temperature for which 
quantum effects become relevant. These are not taken 
into account in the classical Langevin equation How- 
ever, for sufficiently large Josephson junctions a region of 
low temperatures exists for which both thermal and clas- 
sical fluctuations can be neglected on those time scales 
that are experimentally relevant. 

The averaged transport behavior is completely deter- 
mined by the current-voltage characteristic, i.e. the func- 
tional dependence of the averaged voltage on the applied 
dc-strength in the asymptotic limit of large times when 
all transient phenomena have died out. To obtain this 
current-voltage characteristic, we numerically simulated 
10'^ solutions of Eq. ^ from which we estimated the 
stationary dimensionless voltage defined as 



v = m')), 



(3) 



where the brackets denote averages (i) over the initial 
conditions (0(0), <?!>(0), lpq) according to a uniform distri- 
bution on the cube {0(0) € [O,27r],0(O) S [-2,2],(po € 
[0, 27r]}, (ii) over realizations of thermal noise T{t') and 
(iii) a temporal average over one cycle period of the ex- 
ternal ac-driving once the result of the first two averages 
have evolved into a periodic function of time. The sta- 
tionary physical voltage is then expressed as 
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For a vanishing dc-strength, iq = 0, also the average volt- 
age must vanish because under this condition Eq. (2) as 
well as the probability distribution with respect to which 
the average is performed are invariant under the trans- 
formation ((/), ifo) (—0, ipo + I")- For non-zero currents 
io 7^ 0, this symmetry is broken and the averaged volt- 
age can take non-zero values, which typically assume the 
same sign as the bias current io. Apart from this "stan- 
dard" behavior, a Josephson junction may also exhibit 
other more exotic features, such as absolute negative con- 
ductance (ANC) negative differential conductance 
(NDC), negative-valued nonlinear conductance (NNC) 
and reentrant effects into states of negative conductance 
0, §. In mechanical, particle- like motion terms, these 
exotic transport patterns correspond to different forms 
of negative mobility of a Brownian particle. 
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FIG. 1: (color online) Various performance measures of a 
Josephson junction in dependence of the dc-strength are com- 
pared. In panel (a) the stationary voltage v defined in Eq. Q 
is depicted while panels (b), (c) and (d) display the standard 
deviation as a measure the voltage fluctuations, cf. eq. 
((S)), the Peclet number Pe (Eq. ([TTJ) and the two efficiency 
measures rjE (Eq. [7]) and 77s (Eq. |9]), respectively. The re- 
maining system's parameters are: the thermal noise intensity 
D = 0.001, ac driving frequency fl — 0.78, ac driving ampli- 
tude ii — 0.732 and damping constant 7 — 0.143. Three thin 
vertical lines mark the three dc-strengths io ~ 0.0159 (regime 
of absolute negative conductance), io = 0.1114 (Ohmic-like 
regime, chaotic regime, large diffusion, relatively small cur- 
rent fluctuations), io = 0.175 (regime of negative differential 
conductance, very regular motion, small diffusion but rela- 
tively large current fluctuations). For corresponding trajec- 
tories and phase portraits see Figs. [2j [3] and ID 



III. TRANSPORT CHARACTERISTICS 

Apart from the averaged stationary velocity v, which 
presents the basic transport measure, there are other 
quantities that characterize the random deviations of the 
voltage about its average v at large times such as the 
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Here the average is performed with respect to the same 
probability distribution as for v in Eq. ([3]) . This variance 
determines the range 



v{t') e {v ■ 
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of the dimensionless voltage v{t') = 4>{t') in which its 
actual value is typically found. Therefore the voltage 
may assume the opposite sign to the average voltage v if 

ay > V. 

In order to quantify the effectiveness of a device in 
terms of the power output at a given input, several mea- 
sures have been proposed in the literature [H, [H, [13, [H, 
[Tgl . [20j . Here we discuss two of them, which yield consis- 
tent results. For the systems described by Eq. the 
ejjiciency of energy conversion is defined as the ratio of 
the power P = {qv done against an external bias io and 
the input power Pin [21I [22 | , 
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where Pm is the total ac and dc power supplied to the 
system. It is given by pjj 



P,n = 7[v^ + al - D] = -f[{,p^) ^ D]. 



(8) 



This relation follows from an energy balance of the un- 
derlying equation of motion Further, the Stokes ef- 
ficiency is given by the relation 



p. 

^ 1.71. 



D' 



(9) 



where i^^ = jv denotes the Ohmic current, cf. Eq. 
In contrast to the definition of the efficiency of energy 
conversion r]E, the definition of the Stokes efficiency r]s 
does not contain the damping constant 7. We note that 
a decrease of the voltage variance tr^ leads to a smaller 
input power and hence to an increase of the energetic 
efficiency. 

Another quantity that characterizes the effectiveness of 
transport is the effective diffusion coefficient of the phase 
difference (t>(t), describing the spreading of trajectories 
and fluctuations around the average phases. It is defined 
as follows 



= lim 



2t 



(10) 



The coefficient Z?^ can also be introduced via a gener- 
alized Green-Kubo relation [5|]. Intuitively, the diffusion 
coefficient is small and the transport is more effective 
if the stationary voltage is large and the spread of tra- 
jectories is small. The ratio D^/2tt can be considered 
as a velocity characterizing the phase difference diffusion 
over one period. Its relation to the averaged velocity v of 
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FIG. 2: (color online) In the ANC regime at the dc-strength 
Jo = 0.0159 the deterministic trajectories, D = 0, (blue dotted 
lines) for large times (upper panel) approach either a locked 
solution in one of the wells of the periodic potential or a run- 
ning solution that proceeds in the negative direction opposite 
to the positive dc-bias. The remaining parameters are the 
same as in Fig. 1. The sample trajectories of the Langevin 
equation ^ (red solid lines), resulting for the noise strength 
D — 0.001 and the same other parameters, stay close to ei- 
ther of the deterministic trajectories for some time and then 
switch to the neighborhood of another trajectory. The in- 
set of the upper panel displays the locked and the running 
deterministic trajectory and one realization of a stochastic 
trajectory for five periods of the ac driving. The correspond- 
ing phase portraits for the deterministic running solution and 
the stochastic trajectory are shown in lower panels. 



the phase difference determines the dimensionless Peclet 
number Pe defined as 



Pe = 



2^1 W I 



(11) 



A large Peclet number indicates a motion of mainly regu- 
lar nature. If it is small then random or chaotic influences 
dominate the dynamics. 

Fig. 1 depicts the main transport characteristics. 
Panel (a) represents the dependence of the averaged volt- 
age on the dc-strength. It displays ANC for small dc- 




FIG. 3: (color online) In the Ohmic like regime at io = 0.1114 
(other parameters as in the Figs. 1 and 2) the deterministic 
motion for D = (blue dotted line) is chaotic in the limit of 
large times. The phase difference 0(t) increases on average, 
hence leading to normal conductance, see the upper panel. In 
their qualitative appearance the deterministic trajectories are 
not different from the stochastic trajectories (red solid line) 
for the noise strength D — 0.001. The Poincare map of the de- 
terministic system for the same set of parameters is displayed 
in the lower panel. This map represents the phase difference 
modulo 27r and the corresponding velocity of a single trajec- 
tory taken at integer multiples of the ac driving period, see 
e.g. [l^l . It reveals the typical features of a strange attractor. 



strengths and NDC for larger dc-values. Fig. 2ofRef. [8| 
seemingly exhibits a very similar behavior. The present 
ac-strength ii = 0.732 however is somewhat larger than 
the one chosen in Ref. [1]. As a consequence, the present 
set of parameters leads to ANC already for D = 0. The 
deterministic motion then is governed by two coexisting 
attractors, a locked and a non-chaotic running solution, 
see Fig. 2. 

In panel (b) of Fig. 1 the standard deviation cr^ of 
the voltage is depicted as a function of iq. We note 
that its dependence on the dc-strength is rather com- 
plicated and without any immediately obvious relation 
to the averaged voltage of panel (a). Upon closer inspec- 
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FIG. 4: (color online) In the regime of NDC at io = 0.175 
(all other parameters as in the previous Figs.) the deter- 
ministic motion is represented by periodic trajectories (blue 
dotted lines) running in the positive direction. In this pa- 
rameter regime the stochastic trajectories at D = 0.001 (red 
solid line) stay much closer and for much longer times nearby 
a deterministic trajectory than in the cases displayed in Figs 
2 and 3 (see the magnification in the inset and also compare 
with the inset of the upper panel of Fig. 3). This explains 
the small value of the phase diffusion coefficient in this 
parameter regime, see also panel (c) of Fig. 1. Also in the 
phase portrait shown in the lower panel the stochastic tra- 
jectory (red solid curve) remains most of the time in a close 
vicinity of the deterministic trajectory (blue dotted curve). 



tion, one though observes that the vohage fluctuations 
undergo rapid changes yielding a large standard devia- 
tion in the vicinity of zero bias (at small dc-values) and 
in the interval where the regime of the NDC appears. In 
these regimes, the maxima of the standard deviation are 
found, cf. the inset in the upper panel of Fig. 4, where 
one can reveal oscillations in the potential well. Upon 
a further increase of the dc-strength, the voltage fluc- 
tuations decrease because the influence of the periodic 
potential then becomes weaker. Since in this limit os- 
cillations and back-turns of the trajectories become less 
frequent the standard deviation of the voltage decreases. 



In panel (c) of Fig. 1 we demonstrate the influence 
of the dc-strength on the Peclet number which is a non- 
linear function of the bias iq: for small, increasing dc- 
strengths, it first increases and then drops again; the 
rapid growth of the voltage fluctuations around iq ~ 
0.175 is accompanied by an increasing Peclet number. 
In this regime the trajectories of the noisy dynamics stay 
almost always very close to the periodic deterministic at- 
tractor, trajectories bundle closely together, see the up- 
per panel of Fig. 4. 

Finally we study the efficiency of the device. In panel 
(d) of Fig. 1, the two efficiency quantifiers, tje and rjs 
are presented. Both vanish if the averaged voltage is zero. 
In the regime outside of ANC, the energy conversion ef- 
ficiency r;^; is a monotonically increasing function of the 
dc-strength. In the regime of the NDC (in the vicinity of 
io = 0.175, cf. Fig. 1), this efficiency is almost constant 
with the value tje = 0.3. With a further increase of the 
dc-strength it increases and saturates to the value 1 for 
large ig. The Stokes eflaciency r]s attains a local mini- 
mum in the vicinity of the dc-strength ig = 0.18 and is 
always smaller than the efficiency of the energy conver- 
sion r]E- For large io, it also approaches the value 1. In 
the regime of ANC (shown in the inset of panel (d) in 
Fig. 1), both efficiencies are small, of the order lO^'^. In 
Refs. flrl, the rectification efficiency rjn is introduced. 
Adopting this definition to the system we get 



VR = 



-iov + 



-tov + j[{(l)-^) -D] 



(12) 



For the present system, rjn takes both positive and nega- 
tive values. Moreover, both the corresponding efficiency 
of energy conversion rjE and the Stokes efficiency rys when 
evaluated with the input-denominator as given with Eq. 

assume values larger than unity. Therefore, these 
so evaluated three measures 77/j, rjE and rjs are no longer 
suitable to characterize 'efficiency' in the present context 
with an inertial dynamics determined from Eq. 



IV. SUMMARY 

Although Josephson junctions have been studied and 
explored for many years, still new intriguing properties 
are discovered in these particular devices, which also have 
a great potential to impact novel technologies. They be- 
long to the most promising candidates for solid based 
quantum qubits 24]. ANC in Josephson junction has re- 
cently theoretically been predicted ^] and subsequently 
confirmed experimentally [2^. In the present work, we 
continued the study of the main transport characteris- 
tics of such systems. We presented the voltage-current 
characteristic which manifests a regime of ANC and two 
regimes of NDC ( for io « 0.08 and io « 0.15). These 
effects may be realized under various conditions by a 
proper choice of the system's parameters such as tem- 
perature, frequency and ac amplitude and dc strength. 
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Wc revealed that the voltage fluctuations characterized 
by the voltage standard deviation and the phase differ- 
ence diffusion coefficient assume a non-monotonic behav- 
ior as functions of the external load. The voltage stan- 
dard variation exhibits a global maximum in the second 
regime of the NDC, while the diffusion coefficient has a 
global minimum in this regime. Within the ANC regime 
the energetic efficiency is small while in the regime of the 



NDC it takes much larger values. 
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